Noise induced state transitions, intermittency and universality in the noisy Kuramoto-Sivashinsky 

equation 
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Consider the effect of pure additive noise on the long-time dynamics of the noisy Kuramoto-Sivashinsky (KS) 
equation close to the instability onset. When the noise acts only on the first stable mode (highly degenerate), 
the KS solution undergoes several state transitions, including critical on-off intermittency and stabalized states, 
as the noise strength increases. Similar results are obtained with the Burgers equation. Such noise-induced 
transitions are completely characterized through critical exponents, obtaining the same universality class for 
both equations, and rigorously explained using multiscale techniques. 

PACS numbers: 02.50.-r,05.40.-a,05.45.-a,47.54.-r 



Most physical and technological settings are subject to ran- 
dom fluctuations, which are responsible for many intrigu- 
ing and surprising phenomena fV]. These settings are of- 
ten described by model spatially extended systems (SES), 
i.e. infinite-dimensional dynamical systems with space-time 
dependence and some stochastic forcing 1 2] . A widely stud- 
ied example is the transition between different observed sys- 
tem states as the noise strength is increased beyond a critical 
value. For both pure temporal dynamical systems and fully 
nonlinear SES, it is well known that noise-induced transitions 
are due to multiplicative noise, i.e. noise whose amplitude de- 
pends on the fluctuating variable 111 Hi- The presence of an 
additive noise, i.e. noise that does not depend on the state of 
the system, in addition to multiplicative one, has been shown 
to induce other phase transitions L4J while recently it has been 
shown that pure additive noise, i.e. thermal fluctuations, can 
stabilize linearly unstable solutions of SES ISllfifl. However, a 
satisfactory and systematic description of the effects of ther- 
mal fluctuations on SES as well as a quantitative description 
of such effects in terms of critical state transitions is still lack- 
ing. 

In this Letter we report analytical and numerical evidence 
of pure additive noise-induced transitions in SES. As a main 
case study, we consider the noisy KS equation close to the 
primary bifurcation. We observe numerically a number of 
critical transitions by increasing the noise strength, including 
on-off intermittency, a crucial universal feature of many non- 
linear systems close to criticality reflecting a transition from 
order/coherence to a disordered state (hence understanding the 
statistical properties of intermittency is crucial for the charac- 
terization of this transition). Our numerical observations can 
be fully explained in the context of a multiscale theory for 
SES. 

Noise in weakly nonlinear evolution equations.- We con- 
sider the noisy KS equation 



dtu = — (9^ + ^dt)'^ ~ udxU + (j^, 



(1) 



normalized to 27r-domains so that < v = (tt/L)^, where 
2L is the original length of the system, and with either ho- 



mogeneous Dirichlet Boundary Conditions (DBC) or Periodic 
Boundary Bonditions (PBC). Equation ([T]), with and without 
the noise term, has attracted a lot of attention since it ap- 
pears in a wide variety of physical phenomena and applica- 
tions and it also serves as a canonical reference system of SES 
exhibiting spatiotemporal chaos or dissipative turbulence, e.g. 
reaction-diffusion systems and interfacial instabilities in fluid 
flows 10]. 

We shall assume throughout zero-mean solutions, and we 
study a randomly perturbed regime close to criticality by 
slightly increasing the domain size as L = ttVI + so 
that we write u = 1 — and a = ea, where a represents 
the strength of the noise, with e being a bifurcation param- 
eter; if e = all modes, except the neutral one, are sta- 
ble with the system approaching its rest state as t oo, 
and for < e < 1, a bifurcation occurs leading to a fi- 
nite number of linearly unstable modes (in fact [1/ vT^-^J 
of them). The field u can then be projected onto the set of 
eigenfunctions {ek{x)} fork = 1, 2, ... of the linear opera- 
tor C = —d'^ — d^, such that u{x, t) = Uk{t)ek{x), and 
in the limit e ^ 0, only one single mode, namely ui{t), will 
be unstable. We are interested in the dynamics of ui{t) when 
the stable modes, Uk{t) for k >2, are randomly forced, and in 
particular we focus on the case when only the first stable mode 
(U2) is perturbed, so that the noise term in Eq. ([T]) is written as 
^(x, t) = (,2{t)e2{x), where ^2(^) is some uncorrected Gaus- 
sian noise. This "highly degenerate noise" may give rise to 
a stabilization process of the unstable mode iii ||5|]. Typical 
snapshots of the spatio-temporal evolution of Eq. ([T]) subject 
to DBC with cr = 10 and 60 are depicted in Fig.[T] The dy- 
namic evolution of the first mode amplitude, A{t) = \ui{t)\ 
is calculated for different noise strengths and boundary condi- 
tions. 

In the case of DBC, the Probability Density Function (PDF) 
of A{t) calculated for different values of a is shown in Fig.[2l2. 
For a = 15 the amplitude is characterized by finite fluctua- 
tions that never reach zero (state I, top panel in the inset of 
Fig. [2I2), and the PDF can be fitted to a function of the form 
P{A) = A/'A^i exp {-SA^). For a = 35, the noise is strong 
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FIG. 1 : (Color online) Top panels show typical spatio-temporal evo- 
lution of the noisy KS equation solved with DBC for e = 0.1 and 
cr = 10 (a) and 60 (b) at time intervals At = 100 depicted as solid- 
dashed lines. For clarity, the curves are arbitrarily shifted in the ver- 
tical direction. Bottom panels show the corresponding time average 
of t). 
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FIG. 2: (Color online) Numerical results for Eq. ([U integrated on 
a [— 7r,7r] domain with DBC. (a) PDF of the first-mode amplitude 
A(t) = \ui{t) \ for cr = 15 (•), 35 (x), and 51 (green solid line), 
with e = 0.1. Dashed lines correspond to a data fit using P{A) = 
NA^^ exp {—5A^), where the fitted value ai is related to Eq. ([5]). 
The inset depicts typical fluctuations of the amplitude at each of the 
three states discussed in the text, (b) Maximum of the PDF as a 
function of cr^ for different values of e. The solid line corresponds to 
the theoretical solution obtained by solving numerically Eq. 0]) with 
the coefficients of Eq. (|4]). All curves have been normalized to the 
corresponding minimum value at cr = 10. 



enough to alter the behavior of the PDF, shifting its maximum 
position, Amax, which can now reach zero (state II). We char- 
acterize this transition between I and II by computing A^ax 
for different a (cf. Fig.Et), obtaining A^ax ^ - 8501^/^ 
which gives a critical value of ai 29. Finally, as we increase 
cr up to the value a = 51 the first mode is completely sta- 
bilized (state III), and the fluctuations eventually reach zero, 
defining a second critical transition at an ^ 50. Such sta- 
bilization process can be clearly observed when the solution 
is averaged over time (cf. Fig[T]). 



FIG. 3: (Color online) PDF of the waiting times T between two con- 
secutive bursts observed in the on-off intermittent state II correspond- 
ing to the DBC case. We solved both the KS equation, with cr = 35 
and 50, and the Burgers equation, with cr = 9.2 and cr = 7, by using 
e = 0.1. The solid and dashed lines correspond to the numerical 
solutions of the GL model, Eq. ([3]), by using the KS coefficients of 
Eq. ^ with cr = 50, and the corresponding Burgers coefficients with 
cr = 9, respectively. The dotted line is a data fit to P(T) 
with T = 1.50 ±0.01. The inset shows the waiting times PDF in 
the KS equation with cr = 50 and different values for the threshold, 
namely cth = 0.05, 0.1, 0.2, 0.4, and 0.8. 



The middle panel in the inset of Fig. [2I2 demonstrates that 
state II is characterized by an on-off intermittent behavior of 
the amplitude fluctuations. Such intermittency can be charac- 
terized by studying the PDF of the waiting times T between 
two consecutive bursts, defined as large fluctuations above a 
given threshold, i.e. A{t) > Cth 1 9]. Figure [3] shows the nu- 
merical results obtained by using two noise strengths. For 
a = 50, close to the second critical point (an), the PDF of T 
is given by P{T) ~ with r = 3/2. Interestingly, this 
exponent has been ubiquitously found in many other physi- 
cal systems that display avalanche or intermittent dynamics 
close to criticality, including neuronal activity in cortex, elec- 
troconvection of nematic liquid crystals or fluid flow in porous 
media (lit]. For a = 35, far from an, the PDF is exponen- 
tially corrected. These results do not depend on the choice of 
the threshold value Cth (see inset of Fig. O. 

When the system is solved by imposing PEC, the first state 
transition occurs at a :^ 36 [cf. Fig. (H. Interestingly, as 
we increase the noise strength, the second critical transition 
is no longer observed, and the power-law regime of the am- 
plitude PDF increases with an exponent that is asymptoti- 
cally decreasing up to the value around —0.21. Finally, the 
same analysis for both DBC and PEC has been performed 
for the noisy Burgers equation (used, for example, as a pro- 
totype for ID turbulence, albeit without pressure gradient) : 
dtu = {d^ + l)u + e^u + udxU + ea^. We obtain the values 
(Ji 6.4 and an ~ 9.5 for the DBC case, including on-off 
intermittency with the same exponent for the waiting times' 
PDF (see Fig. [3]). In the PBC case, we find ai 9.5, and as 
with the KS equation, the transition to state III is not observed 
either. This critical phenomenon occurring in both models re- 
flects an underlying universal behavior. Our aim now is to 
explain all these numerical results using multiscale analysis 
(singular perturbation theory). 
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□ Numerical result 
— GL model 
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FIG. 4: (Color online) Numerical results for the stochastic KS equa- 
tion ([T]) integrated on a [— tt, tt] domain with PBC. (a) PDF of the 
first-mode amplitude for a = 25 (•), 45 (*), and 65 (+), with 
e = 0.025. Dashed lines correspond to a data fit using P{A) = 
NA^'^ exp [—5'A^). The inset shows the value of the fitted expo- 
nent a2 as a function of compared to the analytical solution (solid 
line) given by Eq.([Tot. (b) Maximum of the PDF as a function of cr^ 
normalized with the value corresponding at a = 10. The solid line 
corresponds to the numerical solution of Eqs. ([8]) and ([9]), by using 
Eq. (H. 



Multiscale theory.- We first analyze the noisy KS equa- 
tion iHqI] . In the limit of e <C 1, the system is close to the 
bifurcation point and Eq. ([T]) has two widely separated time 
scales, corresponding to the (stable) fast and (unstable) slow 
modes. Considering then the behavior of small solutions at 
time scales of 0(e~^), we define u{x^ t) = ev{x, e^t) to trans- 
form Eq. ([T]) to 

dtv = -e-^{dl^dt)v^dtv-e-\d:,v^e-^(je2i2{t), (2) 

where we have assumed highly degenerate noise. 

Eor DBC, the solution can be expanded in the basis 
{ek{x) — Ck sin (kqx)}, where q = tt/L, c^s are normaliza- 
tion constants, and the single dominant mode ui {t) is real and 
belongs to the null space of jC = —d^ — d^. Also, we stipu- 
late that ^(x, t) = /3{t) sin (2qx), with being white noise, 
{I3(t)l3(t')) = 5(t — t').To obtain the dominant mode ampli- 
tude equation we project the field v in Eq. © onto the null 
space of £ to get k;i = VcV, where Vc is the corresponding 
projector to the null space, and onto its orthogonal subspace 
(stable modes) to get w± = {X — Vc)u, where X is the iden- 
tity operator. Equation (|2]), when written for the variables wi 
and w± is of the form of a fast/slow system of stochastic dif- 
ferential equations (SDEs) for which homogenization theory 
applies [8]. By analyzing the corresponding Fokker-Planck 
equation using singular perturbation theory we obtain a closed 
equation for the distribution function of wi from which we 
can read off a one-dimensional stochastic differential equa- 
tion that is valid in the limit e ^ 0. The resulting equation for 
A{t) = \wi \ is given by the GL equation with multiplicative 
Stratonovich noise: 



71 



-1/2688, 72 = 1/48, 73 = 1/24. 



(4) 



Equation (O has been the subject of several studies (e.g. see 
Ref. and references therein), and the corresponding sta- 
tionary PDF for the random variable A is found to be Jsj] 
P{A) = NA^^ exp(-(^A^), with N a normalization con- 
stant, and 

ai(<7) = 2(1+7i<t2)/(73V)-1, S{a) = 72/(73^). (5) 

As noted in Ref. depending on the location of the maxima 
of the above PDF, there may exist different states describing 
the amplitude A. The interesting point is that all the numerical 
states presented before can be achieved by simply changing 
the value of a. First, we observe that as long as ai > the 
maximum of P{A) occurs at a finite value, Amax > 0, and 
then A is characterized by finite fluctuations around a mean 
value (state I). In contrast, for — 1 < ai < 0, the maximum 
is located at zero, A^ax = 0, and the amplitude fluctuates 
intermittently between zero and a finite value (state II). These 
two states are separated by the critical value: 



(71/2-71)-^/^. 



(6) 



Note that for 71 > 0, this transition can only be observed as 
long as 73 > 271, while it is always observed for 71 < 0. By 
using the values of Eq. we therefore obtain ai = 28.4 in 
excellent agreement with the numerical observation shown at 
Fig. [21). In addition, the critical behavior can be characterised 

as Amax = Wi - Cr2|^/2/(criy^) for (J < (Ji, and Amax = 

otherwise, so that Amax and are the order and control pa- 
rameter, respectively, describing the critical transition. By 
solving numerically Eq. ([3]) with the coefficients of Eq. ^ 
for different a, we find very good agreement between analyt- 
ical and numerical results [cf. Fig.[2t]. If 71 < 0, a second 
transition occurs when ai < —1. The PDF cannot be normal- 
ized and it is given by a Dirac delta function, P{A) = S{A), 
describing a completely stabilized state with A = (state III). 
The critical value an for this second transition is: 



(7) 



i = (1 + 71^')^ - 72^^ + ls(TA/3{t), 



(3) 



yielding an = 51.8, in excellent agreement with the numer- 
ical results (cf. Fig. To obtain analytically the statis- 
tical properties of the waiting times T, we assume that in a 
regime close to the critical point (a < an) the initial value 
of A is below a small given threshold Cth, and we ask for 
the probability P{T) that at time T the amplitude reaches the 
threshold for the first time. In this close-to-zero state we can 
neglect the nonlinear term in Eq. ([3]), and we introduce the 
transformation y = log A obtaining y = 1 -\- 7ia^ + j3al3{t) 
withy G (—00, log Cth]. We thus recognise an underlying 
dynamics described by the well-known first-passage prop- 
erties of the random walk (v^ . giving rise in our case to 
the long-time behavior P(T) - T-^/^ exp {-T/Tq) with 
To = [273a/ (1 + 7ia^)]^, from which in the critical point 
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(J = (Jii we recover the numerically observed pure power-law 
(cf. Fig. O. Clearly, the exponent —3/2 will be universally 
observed in any SES whose dominant mode is described by 
Eq. ©. 

Consider now the case with PBC. The solution 
is then expanded in the exponential Fourier basis 
{e/c(x) = c/e exp (i/cgx)}, for k = 0, ±1, dz2, . . . , 
and the single dominant mode has two components: 
ui{t)ei{x) = yi^m.{qx) + zicos(gx). The noise is now 
given as ^(x, t) = pi{t) sin (2qx) + P2{t) cos (2gx), where 
Pi{t) and p2{t) are uncorrected white random variables. By 
applying our multiscale methodology we obtain: 

^1 = (1 + 271^2)^1 - 72^1^2 + 273(7 (8) 

il = (1 + 27i(7")zi - 72^1^2 + 273(7^/32, (9) 

where A{t) = ^y\^ zf, and 71, 72, and 73 are given 
by Eq. (|4]). The stationary joint PDF for the two vari- 
ables, G{yi,zi), is obtained by computing the corresponding 
stationary two-dimensional Fokker-Planck equation, yielding 
G{yi,zi) (X {yf + zf)^'^^'^ exp [—S'{yl + zf)], where a[ and 
6' are obtained from the expressions in Eq. Q by replacing 

71 and 73 with 271 and 273, respectively. To study the be- 
havior of P{A), we move to a polar coordinate system (A, 6) 
and impose the condition G{yij zi)dyidzi = P(A, 0)dAdO, 
getting P{A) ex ^^p (- J'A^), with 

a2{(T) = a; + 1 = (1 + 27ia2)/(273V2). (10) 

We first note that state transitions can only occur iff 71 < 0, 
with the critical values: ai = y^1/2|7i|, and an = [2(|7i| — 
7!)]"^/^. Interestingly, the second transition can only occur 
as long as 73 < |7i|. Otherwise, state III is never observed, 
and the distribution tends to P{A) ~ A^°° as a 00, with 
Qfoo = By using Eq. dU, the first transition occurs 

at (Ji = 36.3, while the second transition cannot be observed 
with a 00 = —0.21, in excellent agreement with the numerical 
results (cf. Fig.|4]). Finally, when this formalism is applied to 
the Burgers equation, we obtain the coefficients 71 = —1/88, 

72 = 1/12, and 73 = 1/6, giving rise to ai — 6.3 and cru = 
9.4 for DBC, and cfi = 9.4 for PBC, in excellent agreement 
with the numerical results. As with the KS equation, we have 

73 > I71 1, and the second transition is not observed for PBC 
either. 

To conclude, we have presented clear evidence of critical 
transitions in SES induced by pure additive noise. We have 
focused on the KS equation and by adding a stochastic forc- 
ing acting on the first stable mode, we have provided a de- 
tailed and systematic investigation of the transitions between 



different states. In particular, by using multiscale analysis for 
SDEs, we have analytically described the different critical- 
state transitions that are undergone by the amplitude of the 
unstable mode, including on-off intermittency and stabilized 
states, that we have also observed numerically in both the KS 
and Burgers equations. Moreover, the critical exponents for 
both SES are the same, and hence they belong to the same 
universality class. This is in accordance with Yakhot's con- 
jecture 1130 . We believe that our results will motivate further 
analytical and numerical studies on the effect of additive noise 
in general SES. 
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